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Abstract 

We propose a method for spectrally filtering radiation in the VUV wavelength range 
by means of a monochromator constituted by a cell filled with a resonantly absorb- 
ing rare gas. Around particular wavelengths, the gas exhibits narrow-bandwidth 
absorbing resonances following the Fano profile. In particular, within the photon 
energy range 60 eV - 65 eV, the correlation index of the Fano profiles for the photo- 
ionization spectra in Helium is equal to unity, meaning that the minimum of the 
cross-section is exactly zero. For sufficiently large column density in the gas cell, the 
spectrum of the incoming radiation will be attenuated by the background cross- 
section of many orders of magnitude, except for those wavelengths close to the 
point where the cross-section is zero. Remarkable advantages of a gas monochro- 
mator based on this principle are simplicity, efficiency and narrow-bandwidth. A 
gas monochromator installed in the experimental hall of a VUV SASE FEL facility 
would enable the delivery of a single-mode VUV laser beam. The design is iden- 
tical to that of already existing gas attenuator systems for VUV or X-ray FELs. We 
present feasibility study and exemplifications for the FLASH facility in the VUV 
regime. 



1 Introduction 

A large portion of experiments using VUV radiation are presently car- 
ried out at synchrotron radiation facilities and VUV SASE FELs. A high- 
resolution VUV monochromator is an important element of such facilities. 
Photo-absorption spectra of Helium are often used for calibration of these 
monochromators. In fact, as is well-known, the transmittance spectrum in 
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the Helium gas exhibits narrow-bandwidth absorbing resonances down to 
a few meV, which follow the Fano profile 

Here we propose a novel method for spectral filtering based on a remark- 
able feature of the Fano interference phenomenon for Helium. In particular, 
within the photon energy range between 60 eV and 65 eV, the correlation 
index of the Fano profile for the photo-ionization spectra in Helium is equal 
to unity, meaning that the minimum of the cross-section is strictly zero [2]. 
We suggest the exploitation of this feature for filtering purposes, using a 
windowless gas-filled cell equipped with differential pumping. In fact, for 
a sufficiently large column density, the spectrum of the incoming radiation 
will be attenuated by the background cross section of many orders of mag- 
nitude, except for those wavelength near the point where the cross-section 
is zero. In this paper we concentrate on physics issues. The design features 
of a gas monochromator system based on this principle should not differ 
too much from the gas attenuator systems at XFELs |13j,|^. 

A remarkable advantage of gas monochromators based on the above prin- 
ciple is simplicity. In fact, such a device would use no optical components, 
and therefore there are no problem with alignment, heat loading and beam 
wavefront perturbations. Moreover, the peak efficiency of a gas monochro- 
mator would be close to 100%. Also, a very narrow bandwidth would be 
granted. Such a device, installed at any VUV SASE FEL facility would enable 
the delivery of single-mode VUV laser beam. In this case, the shapes of the 
single-shot spectra after the monochromator can be identified with the gas- 
cell transmittance calculated from the first principles, and do not depend 
on the detailed intensity distribution of modes of the incoming spectrum. 

We present a feasibility study and exemplifications for the FLASH SASE 
FEL in the VUV regime. The applicability of the gas monochromator setup 
is obviously not restricted to VUV SASE FEL facilities. Synchrotron radiation 
facilities may benefit from this scheme as well. 



2 Principles of the spectral filtering technique based on the use gas cell 

The output radiation pulses from SASE FEL exhibit poor longitudinal co- 
herence. While a single-shot spectrum consists of many spikes, or modes, 
the average spectrum is constituted by a smooth distribution with a FWHM 
typically corresponding to about 1% in the VUV wavelength range. The 
bandwidth of the VUV beam can be thus further reduced by using a gas 
monochromator. A schematic representation of the setup is shown in Fig. 
1. Each spike in the spectrum has an average width of about 0.03% - 0.04% 
for wavelengths around 20 nm, and is thus larger than the bandwidth of 
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Fig. 1. Scheme of the gas monochromator for VUV SASE PEL sources 



the gas monochromator, which allows for a few meV bandwidth. This fact 
enables the delivery of single-mode VUV laser beams. The shape of the 
single shot spectra after the gas monochromator can be described in terms 
of the transmittance of the gas cell, and do not depend on the details of the 
intensity distribution of the spikes in the incoming spectrum. It is therefore 
possible to properly define amplitude and phase for a single-shot spectrum 
after the monochromator and, hence the single-shot radiation pulse shape 
in the time domain. 



The working principle of the gas monochromator is based on the charac- 
teristic of the photo-absorption cross-sections for Helium at energies below 
the second ionization threshold. In general, the interaction with a soft X-ray 
electromagnetic pulse has an ionizing effect for the atom. A first ionization 
threshold, where the remaining electron in the ion is found in the state with 
principal quantum number n = 1, is found at 24.6 eV, while a second ion- 
ization threshold with n = 2 is found at 65.4 eV. Therefore, after interaction 
with a photon with energy between these two thresholds, one finds the 
Helium ion in the ground state Is, and a free electron with kinetic energy 
equal to the photon energy diminished of 24.6 eV. This single open channel 
can be reached either directly or via decay of a doubly excited state, the 
well-known autoionizing state In this case only one channel is open in 
the continuum, meaning that the final state is the Helium ion in the ground 
state. The two paths towards this channels interfere and one is left with 
Fano lines obeying: 
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where o^, = oi,{A) is the background cross-section, q is known as the asym- 
metry index, while the reduced energy S is defined as 



e = H(^i^, (2) 

Er being the energy of the resonance, with a width T. When the energy 
increases, or one considers other gases, new channels open up. For exam- 
ple, above the second ionization threshold for Helium, when the photon 
energy increases between 65.4 eV to 72.9 eV, a second channel is present, 
characterized by the Helium ion with the bound electron with n = 2 (2s 
state), and free electron energy equal to the photon energy diminished by 
65.4 eV. In this photon energy range we have a set of autoionizing reso- 
nances,which interfere with two open channels (Is or 2s). More in general, 
as the photon energy increases below the double ionization threshold at 79 
eV, other channels corresponding to the Helium ion with higher numbers 
of n (nth ionization threshold) can be reached. Then, as explained for ex- 
ample in ll^, transitions to the continuum may or may not interact with the 
discrete autoionization state, and the expression for the cross-section has to 
be modified as: 



iq + £)2 

o = Ob ^ g2 + ' (3) 

where Ob and Oa are the background cross-sections of the interacting and 
non-interacting transitions to the continuum. Defining the total cross-section 
= (7a + Oh, which one would measure without autoionization, Eq. (3) can 
be cast in the form 



where the correlation coefficient is introduced as 



Eq. (4) reduces to Eq. (3) for = 1, corresponding to the Helium resonances 
under study. In this case, there is always one energy value around the reso- 
nance, at £ = -q, where cr = 0. Around that energy value, radiation can pass 
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through the gas cell almost without attenuation. In fact, the photoabsorp- 
tion cross section o is linked to the light attenuation through a gas medium 
of column density HqI, I being the length of the cell and Hq the gas density, 
via the Beer-Lambert law: 



where Iq is the incident intensity, and I{co) is the attenuated intensity of 
the transmitted light at frequency co. By inspecting Eq. (6) it can be seen 
that when the column density is increased radiation is attenuated at all 
frequencies except at 6 = -q. As a result, in the case of the Helium gas, and 
for energies below the second ionization threshold an increase in the gas 
density transforms the gas cell into a bandpass filter. 

2.1 VUV photo absorption spectra of He for photon energies between 60 eV and 



Let us then consider the series of autoionizing resonances converging to the 
n = 2 level of the Helium ion, corresponding 2n + 1 =3 Rydberg series, 
all converging to the second ionization threshold, corresponding to 65.4 eV. 
These three series all begin with the two Helium electrons characterized by 
ni = 2 and n2 = 2, which corresponds to an energy of 60 eV. There are three 
possible choices for the {h, h) quantum numbers: (2s, np) ^Pq, (2p, ns) ^Pq, 
and {2p,nd) ^Pq. However, the levels of the {2s, np) ^Pq and (2p,ns) ^Pq 
wave functions are nearly degenerate, and the doubly-excited states are 
better described in terms of constructive and destructive superpositions of 
these wave functions, finally yielding the three Rydberg series: (sp, 2n+) ^Pq, 
(sp,2n-) ^Pq, and {2p,nd) ^Pq. Since these three series converge (but do not 
reach) the second ionization threshold, all the correspondent doubly excited 
states He* decay to an Helium ion in the ground state. Out of these three 
series, only the first {sp,2n+) will be of interest to us, as the others are too 
weak. 



The energy-dependent background cross-section expressed in Megabam 
(1Mb = lO-^^cm^) is given by \M 



with A the radiation wavelength in Angstrom units. The asymmetry index q, 
the energy of the resonances Er, and the resonance widths T, have been the 
subject of several calculations. The cross-section for the series {sp,2n+) can 
be calculated using Eq. (3) with the parameters given in Table 1, reproduced 
from L2J. 



I{a)) = lo exp[-nolo{co)] , 



(6) 



65 eV 



Ob{A) = -0.05504 - 1.3624 ■ lO'^A + 3.3822 • 10"^ , 



(7) 



6 



Table 1 

Parameters determining the Fano autoionization profiles for the (sp, 2n+) autoion- 



izing series of Helium, following [2]. 



n 


2 


3 


4 


5 


6 


7 


8 


9 


10 


r(eV) 


0.0378 


0.0083 


0.0038 


0.0014 


0.0008 


0.0005 


0.0003 


0.0002 


0.0001 


Eres (eV) 


60.14 


63.655 


64.466 


64.816 


64.999 


65.108 


65.181 


65.229 


65.263 


q 


-2.57 


-2.5 


-2.5 


-2.5 


-2.5 


-2.5 


-2.5 


-2.5 


-2.5 



10- 




Fig. 2. Fano profiles for the {sp^n+) autoionizing series of Helium. The cross-sec- 
tions are calculated following 

2.2 Doppler broadening 



Due to the finite temperature T of the gas, the frequency distribution of 
an ensemble of radiating atoms will experience a spread. Assuming the gas 
ensemble in thermal equilibrium, the Doppler shift is linked to the Maxwell- 
Boltzman probability distribution of velocities of the atoms. As a result, the 
Fano line in Eq. (3) should be convolved with the Doppler-broadening 
profile P [h] 



exp[- 



Mc^jE'-Ef 



!P(£, £'; T) = ' ^ JZL ^ , (8) 

27tlcTE2 



Mc2 
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where k is the Boltzmann constant, T is the temperature, and M is the 
molecular weight, yielding the Doppler-broadened cross-section 



From Eq. (8) we see that the distribution of frequencies turns out to be a 
Gaussian, with a peak at the resonance frequency and a FWHM given by 



Assuming a temperature of 300 K, the Doppler broadening is pressure inde- 
pendent and turns out to amount to about 0.4 meV for He (molecular weight 
M = 4). For the resonances under study, Tu is substantially smaller than the 
natural width F ~ 4 - 10 meV. It follows that the Fano profiles constitute 
a good approximation to the exact profile, which should be obtained by 
convolving the Fano natural line with the Gaussian profile due to Doppler 
shift. However, in order to quantitatively account for deviations of the peak 
efficiency of our monochromatoiO] from 100%, in simulations we account 
for Doppler broadening effects according to Eq. (9). 



3 Amplitude and phase of the gas cell transmittance 

As discussed before, if a monochromatic electromagnetic pulse of intensity 
Jo and frequency co impinges on a cell of length /, filled with a gas with 
density no, the transmitted intensity obeys the Beer-Lambert law, Eq. (6). As 
a result, the modulus of the transmissivity can be defined as 



where o{co) follows a Fano profile convolved with the Doppler-broadentng 
profile, as has been discussed in the previous Section. In this article we 
will assume that the gas used in the monochromator is Helium, and that 
the energy of the EEL photons is within the range between 60 eV and 65 
eV. These requirements ensure that the cross-section has a zero for S = -q, 

^ Doppler broadening is the only effect responsible for the deviation of the peak 
efficiency from 100%. For example, collision broadening is very small. In fact we 
should compare collision times in the order of 100 ns with the typical inverse 
ionization rate in the order of 100 fs. As a result, we expect an influence in the order 
of 0.0001%, which can be neglected. 




(9) 




(10) 



|T| = exp[-noZc7/2] , 



(11) 
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Fig. 3. Schematic of the gas cell 
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Fig. 4. Schematic of the differential pumping system 

meaning that the correlation parameter = 1. In order to obtain a pass- 
band filter, one needs to increase the gas pressure so to absorb radiation 
outside a small bandwidth from the point where & = -q. The pressure in 
the gas cell is limited by capabilities of the differential pumping system, 
and cannot exceed a few Torr in the present design. Fig. 3 and Fig. 4 show 
the schematic of the gas cell, and its integration in the differential pumping 
system. The design presented here is based on |0]-|11], which were used in 
the design of the LCLS gas attenuator. The length of the gas cell is between 
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Fig. 5. Modulus and phase of the transmissivity of Helium around the n = 4 line of 
the (sp, 2n+) ^Pq Rydberg series. The modulus has been calculated according to Eq. 
(11), while the phase is recovered with the help of the Kramers-Kroning relation 
according to Eq. (13). The left plot refers to a column density Mq/ = lO^^cm"-^. The 
right plot refers to a column density HqI = lO^^cm"^. 

3 m and 5 m. Fig. 3, and cannot be easily increased due to space limitations. 
Even with these limitations on length and pressure, the gas cell can provide 
a very large attenuation coefficient. For example, at a cell length of 3 m and 
a pressure of 1 Torr the column density is HqI ~ lO^^cm"^, and the spectrum 
of the incoming radiation is attenuated by the background cross section of 
about 40 dB. 

Let us consider the third (n = 4) resonant line of the (sp, 2n+) ^Pq Rydberg 
series for Helium, and let us calculate its cross-section with the help o|l]Eq. 
(3) and Eq. (7). The modulus of the transmissivity can be found with the help 
of Eq. (11). In order to exemplify the effects of the pressure increase, we select 
two values for the column-density: HqI = lO^^cm"^, and Hq/ = lO^^cm"^. Since 
|r| is known, it is possible to use the Kramers-Kroning relations to recover 
the phase. In fact one can write 



ln[r(a;)] = ln[|r(a;)|] + iO{a)) = -nla/2 + zO(a;) . (12) 

Note that T{a)) = T{-a)) implies that \T{a))\ = \T{-a))\ and that 0(a;) = 
-0(-(i;). Therefore, using Eq. (12) one also has that ln[r((x')]* = ln[r(-(x')]. 
Then, application of Titchmarsh theorem shows that the analyticity of 
ln[|r(Q)|] on the upper complex Q-plane implies that 



oo 

0(a;) = P — -doj , (13) 

n J co'^ - co^ 



^ As already said, in actual simulations we will also account for Doppler broaden- 
ing effects. For simplicity we do not account for them here. 
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Table 2 

Parameters for the nominal pulse mode of operation used in this paper. 



Units 


Undulator period 


mm 


27.3 


K parameter (rms) 




0.89 


Wavelength (fundamental) 


nm 


20 


Charge 


nC 


0.5 


Electron beam energy 


MeV 


560 



A direct use of Eq. (13), with |r| given as in Eq. (11), yields back the phase 
0(0;). By this, one tacitly assumes that ln[T(Q)] is analytical on the upper 
complex Q-plane. This fact, however, is immediately granted by the fact that 
o is proportional to the imaginary part of the refractive index of the medium, 
and it is well known that the refractive index must obey the Kramers- 
Kroning relation Eq. (13). We therefore used Eq. (13) in order to recover 
the phase of the transmittance. More specifically, we took advantage of a 
publicly available Matlab script [jisl] which serves exactly to that end. The 
final result in terms of modulus and phase of the transmissivity T is shown 
in Fig. 5 for the case HqI = lO^^cm"^, left plo10 and HqI = lO^^cm"^, right plot. 
One can see that an increase in the gas density of a factor 10 transform the 
filter in Fig. 5 left, into the few meV bandpass filter in Fig. 5 right. 

The example studied here for the n = 4 resonance case corresponds to a 
pass-band filter with a 2.1 meV FWHM-bandwidth. As said above, such 
bandwidth is narrower than that of a typical EEL mode, and this enables 
the delivery of single-mode VUV laser beams. Note that the knowledge 
of modulus and phase of the transmittance allows for the characterization 
of the single-mode output both in modulus and phase. This is an obvious 
advantage for scientists wishing to model the interaction between single- 
mode VUV pulses and matter. 



4 Simulations of gas monochromator operation at FLASH 

In this section we present a feasibility study of our scheme. We demonstrate 
the potential of our technique using as an example the parameters for the 
FLASH facility in the VUV region, although our setup can be exploited 
by other EEL and synchrotron radiation facilities in the same photon en- 

^ If the reader will compare this plot with Fig. 5 of il2n . the reader will notice 
differences in \T\. This is due to a misprint in Fig. 5 of |12'], where |Tp and not |T| is 
actually plotted. 
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Fig. 6. Electron beam characteristics at the entrance of the setup at the FLASH, 
(upper left) Current profile, (upper right) horizontal and vertical geometrical emit- 
tance. (lower left) energy profile, (lower right) rms energy spread profile. 



ergy range. The main parameters for the FLASH facility operating around 
are listed in Table 2. Full 3D simulations have been performed in order to 
confirm the scheme feasibility. The simulations are performed with the code 
GENESIS 1.3 |14jl , which uses as input the beam parameters obtained in start 
to end simulations (lE]. The beam distributions are shown in Fig. 6 show- 
ing the current profile, the horizontal and vertical normalized emittances, 
the energy profile and the rms energy spread profile. The focusing system 
along the FLASH undulator consists of doublets between each undulator 
segment. The three possible focusing solutions are a doublet structure, a 
FOFO structure and a FODO structure |16l] . In our simulations we imple- 
mented the FODO focusing type. The average betatron function used were 
about 10 m, resulting in the beam rms size shown in Fig. 7 as a function of 
the distance inside the undulator. 



In the following we will consider different resonances, and different possible 
settings of FLASH, in order to analyze several configurations of the gas 
monochromator setup. 
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Fig. 7. Electron beam horizontal (left plot) and vertical (right plot) rms sizes as a 
function of the distance inside the undulator. 
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Fig. 8. Left plot: the SASE spectrum at saturation overlaps the first resonance. Right 
plot: the filter bandwidth for a column density of no/ = lO^^cm"^. 



4.1 First resonance at 60.14 eV 



We begin considering the first autoionization profile n = 2 for the (sp, 2n+) 
autoionizing series of Helium. FLASH is tuned at saturation, that is we 
consider three active undulator modules corresponding to a magnetic length 
of 15 m, so that the SASE average bandwidth overlaps the first resonance, 
as shown in Fig. 8 left. The right plot in Fig. 8 shows, instead, the bandpass- 
filter profile in terms of \T\^ for a column density of nol = lO^^cm"^. 

As one can see, in this case the filter bandwidth is around 23 meV FWHM, 
which is narrower, but in the same order of magnitude of the bandwidth 
of a single SASE spike. As a result, we should expect a dependency of the 
output spectra on the input intensity, which changes from shot to shot. The 
envelope of the output spectra will follow, however, the distribution in Fig. 
8, right. 

The input spectra and power are shown respectively in the left and right 
plots of Fig. 9. 
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Fig. 9. First resonance mode of operation. Pulse spectrum (left plot) and power 
(right plot) before the gas cell. Grey lines refer to single shot realizations, the black 
line refers to an average over one hundred realizations. 
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Fig. 10. First resonance mode of operation. Pulse spectrum (left plot) and power 
(right plot) after the gas cell. Grey lines refer to single shot realizations, the black 
line refers to an average over one himdred realizations. 

The effect of the filter is shown in Fig. 10, where we plot the output spectra 
(left plot) and powers (right plot). The shot-to-shot fluctuations in the pro- 
files of both intensity and spectrum are a consequence of the fact that the 
filter bandwidth is of the same order of a single spike profile. 



It should be noted that the theoretical calculations used to calculate the 
overall transmission function of the gas cell are based on the validity of the 
perturbation theory approach. In our case the transverse size of the photon 
beam at the cell is about 1 mm (thus fitting into the orifices aperture of 3 mm 
diameter. Fig. 3). The energy in the VUV pulse can be evaluated from Fig. 9, 
and is about 0.3 mj. At VUV wavelengths, this means that a pulse contains 
about 10^^ photons. The photon density is thus lO^^cm"^. It should be noted 
that the resonance line is about two orders of magnitude narrower(10"^) 
compared with the SASE bandwidth (5 ■ 10"'^), while the maximum cross 
section reaches about 10 Mb. As a result, ionization takes place mainly due 
to the background cross section, of about 1.5 Mb. Recalling that 1 Mb is 
lO^^cm^ one can confirm the validity of the perturbation theory approach 
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Fig. 11. Left plot: the SASE spectrum at saturation overlaps the second resonance. 
Right plot: the filter bandwidth for a column density of UqI = lO^^cm"-^. 

used to calculate the overall transmission function of the gas cell. 

This reasoning is valid for a single radiation pulse. Limitations for time 
diagram of facility operation can be estimated with the help of flTt], where 
space-charge effects are studied, and an estimation of ambipolar diffusion is 
presented. Based on this study we can estimate that the single bunch mode 
operation with a 10 Hz repetition rate at FLASH will not be problematic, 
while for other facilities it should be possible to stretch this range up to 
about 100 Hz. However, considering, for example, the FLASH multi bunch 
mode of operation (1 MHz rep rate within 1 ms macro pulse, with 10 Hz 
macro pulse repetition rate) further studies on the diffusion process in the 
gas cell are needed, which are left for future work. 



4.2 Second resonance at 63.655 eV 



We now consider the autoionization profile n = 3 for the (sp, 2n+) autoioniz- 
ing series of Helium. FLASH is tuned at saturation, so that the SASE average 
bandwidth overlaps the second resonance, as shown in Fig. 11 left. The right 
plot in Fig. 11 shows, instead, the bandwidth of the filter in terms of |Tp for 
a column density of nol = lO^^cm"^. 

The input spectra and power are shown respectively in the left and right 
plots of Fig. 12. 

The effect of the filter is shown in Fig. 13, where we plot the output spectra 
(left plot) and power (right plot). The shot-to-shot fluctuations in the profiles 
of both intensity and spectrum are almost absent (although, for a single 
mode, we have a power fluctuation near to 100%). In fact, in this case the 
filter bandwidth is around 5.3 meV IWHM, which is already much narrower 
than the bandwidth of a single SASE spike. As a result, we are approaching 
the single-mode case. 
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Fig. 12. Second resonance mode of operation. Pulse spectrum (left plot) and power 
(right plot) before the gas cell. Grey lines refer to single shot realizations, the black 
line refers to an average over one hundred realizations. 
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Fig. 13. Second resonance mode of operation. Pulse spectrum (left plot) and power 
(right plot) after the gas cell. Grey lines refer to single shot realizations, the black 
line refers to an average over one hundred realizations. 



4.3 Third resonance at 64.466 eV 



We now turn to consider the third autoionization profile for n = 4 for 
the {sp, 2n+) autoionizing series of Helium. FLASH is tuned at saturation, 
so that the SASE average bandwidth overlaps the third resonance. In this 
paragraph we focus on the transition between bandstop to passband filter, 
as the gas column density increases. 

As discussed before, when the gas pressure is decreased, the filter ceases 
to behave as a bandstop filter. This feature is shown by the comparison of 
Fig. 14, Fig. 15 and Fig. 16, which show, respectively, the output spectra and 
power for a column density of UqI = lO^^cm"^, nol = 5 • lO^^cm"'^ and nol = 
lO^^cm"'^. In Fig. 17 we show the average power spectra for a column density 
varying from nol = lO^^cm"^ to nol = lO^^cm"^ in steps of nol = lO^^cm"'^. 

In the case under study, the resonance width is only 2.1 meV FWHM. For 
small gas colunm densities, the filter behaves as a bandstop filter as can 
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Fig. 14. Third resonance mode of operation. Pulse spectrum (left plot) and power 
(right plot) after the gas cell. Grey lines refer to single shot realizations, the black 
line refers to an average over one hundred realizations. Plots refer to a column 
density of HqI = 10^^ cm"-^. 
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Fig. 15. Third resonance mode of operation. Pulse spectrum (left plot) and power 
(right plot) after the gas cell. Grey lines refer to single shot realizations, the black 
line refers to an average over one hundred realizations. Plots refer to a column 
density of UqI 



5 • lOi^cm-2. 



be seen from Fig. 17. The inset of Fig. 17, or the left plot in Fig. 16 shows 
the output narrow bandwidth. As the gas pressure increases, we have an 
attenuation of the spectrum of many orders of magnitude, but the peak in 
the spectrum remains the same for all pressures. In other words, the peak 
efficiency of our monochromator is close to 100%, the difference being only 
due to Doppler broadening, which, as we demonstrated before, constitutes 
a very small effect. For comparison, a conventional monochromator exhibits 
a peak efficiency of about 5% - 10%. As for the case of n = 3, the shot-to-shot 
fluctuations in the profiles of both intensity and spectrum are almost absent 
(although, for a single mode, we have a power fluctuation near to 100%), 
and we basically deal with a single-mode laser pulse. 

Except for the amplitude, the characteristics of the single-mode pulse are 
fixed by the filter transmittance and not by the SASE spectra, which changes 
from shot to shot. This can be shown in the frequency domain by comparing. 
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Fig. 16. Third resonance mode of operation. Pulse spectrum (left plot) and power 
(right plot) after the gas cell. Grey lines refer to single shot realizations, the black 
line refers to an average over one hundred realizations. Plots refer to a column 
density of UqI = 10^^ cm"-^. 
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Fig. 17. Third resonance mode of operation. Pulse average spectra for a column 
density varying from tiqI = lO^^cm"-^ to UqI = lO^^cm"-^ in steps of fiol = lO^^cm"^. 
The lines refers to averages over one hundred realizations. The inset shows an 
enlargement of the bandpass range. 

for example, the spectra in Fig. 16 with the square modulus of the right plot 
in Fig. 5. In particular, we compared the average spectrum in Fig. 16, left, 
with the square modulus of the transmittance calculated with Eq. (11) and 
Eq. (13) (and shown in Fig. 5). The result is plotted in Fig. 18, left, where all 
plots are normalized to unity. It is also interesting to make such comparison 
in the time domain. This can be accomplished by comparing the right plot 
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Fig. 18. Third resonance mode of operation. Left plot: comparison between the 
average spectrum in Fig. 16, left, with the square modulus of the transmittance 
calculated with Eq. (11) and Eq. (13). Right plot: comparison between the average 
power profile in Fig. 16, right, with the square modulus of the Fourier transform of 
the transmittance calculated with Eq. (11) and Eq. (13). All plots are normalized to 
unity. The gas column density is tiqI — lO^^cm"-^. 

in Fig. 16 with the Fourier transform of the transmittance. In particular 
we compared the average power profile in Fig. 16, right, with the square 
modulus of the Fourier transform of the transmittance calculated with Eq. 
(11) and Eq. (13). The result is plotted in Fig. 18, right, where all plots are 
normalized to unity. 

4.4 Production of doublet spectral lines corresponding to the second and to the 
third resonance 

Up to now, we considered cases where FLASH was tuned at saturation, 
so that the SASE average bandwidth overlapped with one resonance at a 
time, n = 2, 3 or 4. However beyond saturation, the bandwidth of the SASE 
spectrum continues to increase, due to a deterioration of the coherence 
time. It is therefore possible to set the operating point of FLASH well after 
saturation, that is after six active undulator modules, corresponding to a 
magnetic length of 30 m. In this way the SASE spectrum can be made 
overlap two (or, in principle, more) resonances, yielding a doublet (or, in 
principle, a multiplet) structure. The overlap between the spectrum of an 
incoming radiation pulse (single shot) and the relevant resonances for n = 3 
(second) and n = 4 (third) are shown of Fig. 19. 

The input spectra and power are shown respectively in the left and right 
plots of Fig. 20. 

The output spectra and power are shown respectively in the left and right 
plots in Fig. 21. The effect of the superposition with the double resonance 
is clear from the analysis of the left plot. In is interesting to note that since 
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Fig. 19. The SASE spectrum after saturation overlaps the second and the third 
resonance. Here a single shot spectrum after six undulator modules is shown 
together with the Fano profiles relevant to the two resonances involved. 
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Fig. 20. Pulse spectrum (left plot) and power (right plot) before the gas cell. Grey 
lines refer to single shot realizations, the black line refers to an average over one 
hundred realizations. 

there are now two frequency filters centered at different central frequencies, 
a beating between the two waves at slightly different frequencies takes 
place, yielding a power distribution modulated in the infrared range at a 
few microns wavelength. This effect is evident from the analysis of the inset 
of Fig. 22, where we show the output power form a single shot pulse. This 
effect is similar to that considered in [IS*]. In that paper we were considering 
the possibility of self-seeding at two or more frequencies in the hard X-ray 
region, and the beating between the two waves was also present in that case. 
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Fig. 21. Pulse spectrum (left plot) and power (right plot) after the gas cell. Grey 
lines refer to single shot realizations, the black line refers to an average over one 
hundred realizations. 
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Fig. 22. Single-shot power distribution. The two frequency filters centered at dif- 
ferent central frequencies are responsible for a beating between the two waves at 
slightly different frequencies, yielding (see inset) a power distribution modulated 
in the infrared range at a few microns wavelength. 

5 Conclusions 

We propose a new method to obtain a spectral filter in the VUV wavelength 
range using the Fano interference phenomenon for the VUV photo absorp- 
tion spectra of Helium. Our gas monochromator combines a very narrow 
bandwidth, down to a few meV and a very high peak efficiency of 99% with 
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a much-needed experimental simplicity. No optical elements nor alignment 
are required. The applicability of our scheme is limited to the wavelength 
range around 20 nm. In this range, however it is possible to completely de- 
fine the output VUV pulse characteristics, both in modulus and phase. This 
opens up the possibility of full characterization of the VUV radiation pulse 
in the time domain, which constitutes a powerful advantage for materials 
scientists wishing to model radiation-matter interactions. 
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